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t-H Triply heavy baryons are very interesting systems analogous to heavy quarkonia, but are difficult 

to access experimentally. Lattice QCD can provide precise predictions for these systems, which 
can be compared to other theoretical approaches. In this work, the spectrum of excited states of 
the £lbbb baryon is calculated using lattice NRQCD for the b quarks, and using a domain-wall 
action for the u, d and s sea quarks. The calculations are done for multiple values of the sea-quark 
masses, and for two different lattice spacings. The energies of states with angular momentum 
up to J = 7/2 are calculated, and the effects of rotational symmetry breaking by the lattice are 
analyzed. Precise results are obtained even for the small spin-dependent energy splittings, and the 
contributions of individual NRQCD interactions to these energy splittings are studied. The results 
are compared to potential-model calculations. 
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1. Introduction 

The bbb system can be viewed as the baryonic analogue of the bottomonium system. Like 
bottomonia, bbb baryons are governed by multiple well-separated energy scales and are therefore 
amenable to the description with effective field theories [1]. Baryons exhibit the SU(3) gauge 
symmetry of QCD more directly than mesons, and are sensitive to the resulting genuine three -body 
forces. Triply heavy baryons probe the three-quark interactions at relatively short distances, where 
contact with perturbation theory can be made [2]. 

The aim of the work reported here is to complement the perturbative QCD calculations and 
other theoretical studies of triply heavy baryons with nonperturbative lattice QCD calculations. 
In Ref. [3], the mass of the ground-state Q-bbb baryon was calculated using lattice QCD to be 
14.371 ± 0.004 stat ± 0.01 1 syst ± 0.001 exp GeV. This was followed by a lattice QCD calculation of 
Q.bbb excited states in Ref. [4], which is summarized here. Like the bottomonium spectrum, the bbb 
spectrum features a hierarchy of radial/orbital excitations and spin-dependent energy splittings. In 
Ref. [4], the energies of ten bbb excited states where calculated with high precision, resolving even 
the smallest spin-dependent energy splittings. The calculation includes dynamical u, d, and s sea 
quarks, and was done at two different lattice spacings of approximately 0.11 fm and 0.08 fm. For 
lattice spacings of this order, the b quarks can be implemented accurately using lattice NRQCD 
[5]. With lattice NRQCD it is also possible to study individually the effects of the different spin- 
dependent interactions in the effective field theory on the bbb energy splittings. 

2. Lattice actions 

The lattice gauge-field configurations used in this work to perform the path integral were 
generated by the RBC/UKQCD collaboration and are described in Ref. [6]. These configurations 
include the vacuum-polarization effects of the light and strange quarks. The Iwasaki discretization 
was chosen for the gauge action, and the u, d, and s quarks were implemented using a domain-wall 
action that preserves chiral symmetry even at nonzero lattice spacing. The main parameters of the 
ensembles are given in Table 1 . 
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Table 1: The lattice sizes, the gauge couplings (/3 = 6/g 2 ), the quark masses, and the corresponding lattice 
spacings and pion masses. 
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In the following, we denote the ft-quark field by y/\ Similarly to Ref. [5], the lattice NRQCD 
action is written as 

wiy cv) - (i - °f ) (> - f (' - f )"(■ - 

(2.1) 

where t/4 are the temporal gauge links, and Hq and 8H are given by 

A( 2 ) 

#0 — —Z , 

2mt, 

8H = - Cl ^i + C2-^fV-E-E-V)-c 3 / T (T-fVxE-ExV)- C 4^ff-B 

8ml smf V / 8m£ V / liny 

+ c 5 ~Wa c 6 T~ 

24mb I6nmf 

-CT-^fA^.a-Bl-cg^r/A^, CT-fVxE-ExVU-c 9 ^<T-(ExE). 

(2.2) 

The different terms in the NRQCD action are suppressed by different powers of v, the average 
speed of the b quarks inside the hadron. The leading term, Ho, is of order v 2 , and gives the dom- 
inant contribution to radial/orbital energy splittings in the bb and bbb systems. The terms with 
coefficients ci 2,3,4 are the relativistic corrections of order v 4 , and the terms with coefficients £7,8,9 
are relativistic corrections of order v 6 (only the spin-dependent order-v 6 terms are included). The 
terms with coefficients C5 6, which contain powers of the lattice spacing, reduce discretization er- 
rors and are not present in the continuum NRQCD action. Because spin splittings first arise at 
order v 4 through the operators ff-|VxE-ExVj and a ■ B, the matching coefficients C3 and C4 
of these operators were tuned nonperturbatively to achieve high precision [4]. The other matching 
coefficients were set to their tree-level values, c/ = 1. 

3. Construction of bbb operators 

The bbb energies can be extracted from Euclidean two-point functions of interpolating oper- 
ators Q.^ with the desired quantum numbers of the bbb states on the lattice. These operators were 
constructed starting from Gaussian-smeared b quark fields 



Waa 



2n s ) 



(3.1) 



where a = 1,2,3 is the color index, a =t,4- * s tne s P m index, A^ 2 ^ is a gauge-covariant lattice 
Laplace operator, and the smearing radius was chosen to be rs ~ 0.14 fm. To obtain a nontrivial 
spatial structure, up to 2 gauge-covariant derivatives were then applied to obtain the following 13 
quark building blocks, 

Waal = Waa, 

Waal = (V. v W)aa, Waa3 = (Vy \j/)aa, WaaA = (V z W)aa, 

Waa5 = {VxVxW)aa, Waa6 = (Vy V* W) aa , ... , Waal3 = (V zW)aa- (3.2) 
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Using Clebsch-Gordan coefficients, these building blocks were then combined to baryon operators 
Q. J m that would have a definite total angular momentum J in continuous space: 

^m= {JM L ^ m L\S,m s } r ijka p r (L,m L ,S,m s ) s a bc %ai¥bj5j¥crk- (3.3) 
m L , mi- 
There are multiple ways of combining the spin indices a,j6,y and the derivative indices i,j,k to 
a given total spin S and total orbital angular momentum L, leading to different types of operators 
classified by permutation symmetries [7], With two derivatives, the largest possible value of 7 is Z. 

Next, in order to account for the breaking of rotational symmetry by the lattice, linear com- 
binations of the different m-components of the operators Q J m were formed, such that the resulting 
operators transform in irreducible representations, A, of the double cover of the octahedral group 
2 0: 

# = (3-4) 

m 

The subduction coefficients S^/" can be found in Ref. [7]. The index r = l...dim(A) denotes the 
row of the irrep. The full procedure described in this section produces seven different bbb operators 
in the H g irrep, three different operators each in the G\ g and Gig irreps, and one operator each in 
the H u and G\ u irreps [4], where the subscripts g/u denote even/odd parity. 



4. Data analysis 

Using the operators from Eq. (3.4), the row-averaged two-point functions 

i dim(A) 



l 



were computed. Here, the indices i,j label the different operators constructed in a given irrep A. 
For example, in the G\ g irrep, there are three different operators which were constructed using the 
following values of (L, S, J): (0, |, |), (2,|,^), and (2, |, |). The numerical results for (4.1) were 
then fitted using the functions 



Q$W = lAiX^\ (4.2) 

n=\ 



where the number of exponentials, N, was chosen equal to the number of operators for each irrep. 
The full spectral decomposition of Cfj{t) also contains an infinite tower of higher excited states, 
but the starting value of t was chosen large enough such that the contributions from these states are 
negligible. 

Sample fit results for E„' and A$ are shown in Fig. 1. The value of J from which each 
operator £2 A W was subduced is also shown. It turns out that for each energy level E ( n A) , the am- 
plitudes A^f are large only for one value of /, and this value of J can therefore be assigned to the 
energy level E„ . The operators retain a strong "memory" of the original value of / [7], and also 
L and S [4]. Approximately matching energy levels with J = | are seen to appear in both the H 
and G2 irreps, while approximately matching energy levels with / = Z appear in all three irreps, as 
expected. 
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Figure 1: 



L.O 

0.8451(17) 0.5 
0.0 

1.0 

0.8290(17) 0.5 
0.0 



...D 



0.8257(18 



G 2 o 



L.O 
0.5 
0.0 



i= 1 
r- 5 



Fitted energies aE„ and relative amplitudes A^f /A^} in the H g 



Gig, and G2 g irreps. The 



values of J from which the operators £2 A W were subduced are indicated below the plots and by the colors. 
The data shown here are for a « 0.08 fm, am u j = 0.004. 
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Figure 2: Examples of chiral extrapolations: E\(l 
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It is shown in Ref . [4] that the small splittings of matching energy levels between the different ir- 
reps become smaller when the lattice spacing is reduced, demonstrating the restoration of rotational 
symmetry. To get the best estimates of the continuum energy levels, simultaneous fits were per- 
formed across the different irreps, using common energy parameters for the previously identified 
matching levels [4]. With the values of J assigned, the energies can then be relabeled as E n {J p ) by 
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Figure 3: Final results for the bbb energy splittings with respect to ground state ), compared to the 

potential-model calculation of Ref. [8]. For the lattice QCD results, the inner shaded bands give the statis- 
tical/fitting/scale setting uncertainty, while the outer shaded bands also include the systematic uncertainty. 
Because the results for the different energy levels are highly correlated, splittings between nearby energy 
levels can be computed with much smaller absolute uncertainties (see Fig. 4). 
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Figure 4: Final results for the bbb energy splittings with respect Xa E\(^ ) (left) and E\{^ ) (right). The 
inner shaded bands give the statistical/fitting/scale setting uncertainty, while the outer shaded bands also 
include the systematic uncertainty. 



increasing value for each J p channel. Note that, due to the use of NRQCD, only energy differences 
are physical. The energy differences computed on the various gauge-field ensembles were then 
extrapolated to the physical value of the pion mass. These extrapolations were done simultaneously 
for the two different lattice spacings, as shown in Fig. 2. The results at the finer lattice spacing and 
at the physical pion mass can be taken as the final values for the physical bbb spectrum. The 
remaining systematic uncertainties were estimated individually for each energy splitting in Ref. [4] 
by considering the contributions of the various NRQCD interaction terms. 
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In Fig. 3, the final results for the bbb energy splittings with respect to the ground state are 
compared to the potential-model calculation of Ref . [8] , where the potential is a sum of two-body 
interactions and includes spin-spin interactions, but does not include spin-orbit or tensor interac- 
tions. For the gross structure of the spectrum, the agreement of this model with the lattice QCD 
results is remarkably good. However, because the model has no spin-orbit or tensor interactions, 
it predicts E 2 {\ + ) = £ 3 (| + ) = £i(§ + ) = E x {\ + ), E A {\ + ) = £ 2 (§ + ), and E x (f) = As 
can be seen in Fig. 4, all of these degeneracies are lifted in full QCD. 



5. Contributions of individual NRQCD interactions 

An interesting question is the following: what are the effects of the different spin-dependent 
operators in the NRQCD action, Eq. (2.2), on the bbb energy splittings? This was investigated in 
Ref. [4] on one ensemble of gauge configurations (a « 0.1 1 fm, am u j = 0.005) by setting different 
subsets of the NRQCD coefficients, c ; , to zero, and recalculating the bbb spectrum in each case. 
The results are shown here in Fig. 5. In case (a), all spin-dependent interactions are turned off. Up 
to small corrections caused by the breaking of rotational symmetry on the lattice, this leads to the 
degeneracies E 2 {\ + ) = £ 3 (| + ) =E l {\ + ) =E,{\ + ), E A {\ + ) = E 2 {\ + ), and E x {\~) = 
which are also present in the potential model of Ref. [8]. 

In case (b), the operator a ■ B is turned on. Interestingly, the biggest effect is seen for negative- 
parity energy levels, where this operator introduces a splitting of £1(5 ) — E\{\ ) = —12.97(45) 
MeV. Note that potential models with both spin-spin and tensor interactions, but without spin-orbit 
interactions, predict £1(5 ) — E\ ( | ) = for Q. baryons [9] , which suggests that the splitting seen 
here is of spin-orbit type. 
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Figure 5: Dependence of the energy splittings E„(J + ) — E\(l ) (left) and E\{\ )—Ei{\ ) (right) on the 
coefficients of the spin-dependent NRQCD interactions [see Eq. (2.2)]. For each energy level, five different 
choices of coefficients are shown. These results are from the ensemble with a ~ 0.1 1 fm, am u 4 = 0.005. 
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As shown in case (c), the operator a ■ (V x E — E x V) has a large effect on both the positive- 
parity and negative-parity bbb energy splittings, which is of opposite sign to the effect of a ■ B. 
For the multiplet {E2(h + ),E^(i + ),Ei(j ),Ei(l )}, which in quark models has L = 2, 5 = |, 
and a totally symmetric spatial wave function, the energy shifts due to G ■ (V x E — E x V) are 
approximately proportional to 2L-S = /(/+l) — L(L + 1) — 5(5+1), which is characteristic for 
a spin-orbit interaction [10]. 

In case (d), both a ■ B and (T-(VxE-ExV) are turned on, and the resulting shifts of the bbb 
energy levels are approximately equal to the sum of the shifts from (b) and (c), but some deviations 
from linearity can be resolved [4]. Finally, in case (e), the spin-dependent order-v 6 interactions 
are also included, restoring the full action as used in the main calculations of this work. The bbb 
spin splittings change by up to 30% when these terms are included. A similarly large effect of the 
order-v 6 terms has also been observed for the bottomonium spin splittings [1 1]. 

6. Conclusions 

The bbb system is an excellent benchmark for nonrelativistic effective field theories and po- 
tential models. Here, a precise lattice QCD calculation of bbb excited states up to J = | was pre- 
sented. Spin-dependent energy splittings that were neglected in potential models were resolved in 
this work, and it was shown what contributions these splittings receive from the different NRQCD 
interactions. The author hopes that these results will stimulate further work on triply heavy baryons 
using a variety of approaches. 
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